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Abstract. We have obtained deep color-magnitude diagrams (CMDs) of the Sagittarius dwarf irregular galaxy with the 
Advanced Camera for Surveys (ACS) on board the Hubble Space Telescope. The new diagrams show for the first time the 
star formation history of SagDIG. The young main-sequence is traced down to ;?if606w - 27.5, well separated from the older He 
burning blue loop stars. The wide color range spanned by the blue and red super-giants is consistent with Padua evolutionary 
models with Z ~ 0.0004. A well-defined red giant branch indicates that the galaxy is dominated by a stellar population older 
than 1 Gyr, for which we have estimated a metallicity in the range [Fe/H] = -2.2 to —1.9 depending on the adopted redden- 
ing. The identification of several C stars on a luminous asymptotic giant branch, as well as the presence of a conspicuous red 
clump with an elongated structure, are indicative of an extended star formation epoch between 1 and 10 Gyr ago. Perhaps most 
importantly, we have identified a small population of stars on a genuinely old red horizontal branch at m F6 06w= 25.70 ± 0.05, 
which sets the first epoch of star formation in SagDIG back to ~ 9 - 10 Gyr ago. This shows that SagDIG has been able to 
maintain a very low mean metallicity over a 10 Gyr life time, and that all Local Group dwarf galaxies share an ancient (> 10 
Gyr) initial episode of star formation. Thus, at the moment, IZw 18 remains the only exception to this general trend, favoring 
the so-called "young galaxy" scenario where a very low metallicity galaxy is not necessarily old. However, a re-analysis of 
ACS observations of IZw 18 suggests that an older RGB population may be present also in IZw 18. 

In order to address the interplay between stars and Inter-Stellar Medium, we compare the distribution of atomic hydrogen 
with that of stellar populations of various ages. The youngest stars are found in clumps located: (i) near to, but not coincident 
with, the highest column densities of H i, and (ii) in tails reminiscent of propagating star formation. We find, however, little 
correlation between stellar populations older than ~ 1 Gyr and the H i. In particular we fail to find any generation of stars that 
preferentially lie within the large H i hole. We therefore exclude the possibility that the ring-like H i distribution is the result 
of multiple supernova events. Alternative scenarios have to be explored for the ring-like H i morphology of SagDIG, which 
amongst dwarf irregular galaxies, is not unusual. 

Key words, galaxies: individual (SagDIG) — individual (IZw 18) — galaxies: structure — galaxies: dwarf — galaxies: pho- 
tometry — galaxies: Local Group — galaxies: fundamental parameters — galaxies: ISM — Galaxy: stellar content — stars: 
formation 



1. Introduction 

The Hubble Space Telescope has made it possible to resolve 
the stellar populations of dwarf irregular galaxies (dlrr) into in- 
dividual stars both in the Local Group (e.g. Dohm-Palmer et al. 
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nW71[rW8ll2"nfJ2l Tolstoy et al. lT9951 Dolphin 120001 Dolphin 
et al. I200T1 120021 120051 Skillman et al. 120051 and beyond (as 
in the case of the blue-compact dwarf IZw 18 at a distance be- 
tween 10-15 Mpc; Aloisi et al. [l999l Izotov & Thuan l2004l . In 
all cases, color-magnitude diagrams (CMD) can be constructed 
to estimate their distance and metallicity and constrain the re- 
cent star formation history (SFH) of the systems. In the nearest 
galaxies, a wealth of evolutionary details can be observed, oth- 
erwise hampered by poor resolution in ground-based studies. 
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The quest for very metal-poor galaxies among dwarf ir- 
regulars started soon after the Searle and Sargent ( 119721 dis- 
covery that the interstellar medium (ISM) of I Zw 18 has an 
oxygen abundance that is ~ 1/50 of the solar one. As Searle 
& Sargent pointed out, this is compatible with two scenarios: 
either the galaxy is now experiencing its first star-formation 
(SF) episode, or, if it contains ~ 10 Gyr old stars, its SF must 
have proceeded in strong bursts separated by quiescent phases. 
The first hypothesis (a young galaxy) is relevant in the context 
of cosmological models that assume hierarchical formation of 
structures in a cold dark matter dominated Universe. In that 
scenario, dwarf dark matter halos (~ 10 8 M ) would be the first 
objects to condense out of the Hubble flow, and it is therefore 
important to understand whether they already contained stars 
before merging into larger galaxies. In fact, IZw 18 is a diffi- 
cult object to study because of its distance. We may ask our- 
selves an equivalent, possibly easier to answer, question: are 
there any galaxies with a very metal-poor young population 
and, at the same time, a ~ 10 Gyr old generation of stars ? The 
existence of such old, metal-poor stellar systems would imply 
that very low metal abundances are not inconsistent with an an- 
cient epoch of galaxy formation. In order to answer this ques- 
tion, we have undertaken an HST/ACS study the Sagittarius 
dwarf irregular galaxy (also known as SagDIG or UKS 1927- 
177). The interest of SagDIG arises from the presence of on- 
going star formation in an extremely metal-poor environment. 
The fact that SagDIG is a nearby (^1.1 Mpc) galaxy makes it 
an ideal target to test locally the young galaxy hypothesis. 

The most recent photometric studies of SagDIG are those 
of Karachentsev et al. (119991 . Lee & Kim ( T2000> . and Momany 
et al. 120021 hereafter M02). Lee & Kim and Karachentsev and 
coll. concluded that SagDIG has the lowest stellar metallicity 
(up to [Fe/H] 2.45) among Local Group star-forming galax- 
ies. M02 presented deep BVI photometry of SagDIG obtained 
from ESO NTT data (we refer the reader to that paper for a 
full summary of previous work). In M02 a differential redden- 
ing scenario was suggested to explain the relatively red color 
of main-sequence stars. Our hypothesis was that some inter- 
nal reddening (in addition to the foreground reddening) affects 
the young stars in star-forming regions. The revised metallicity, 
[Fe/H] = -2.1 + 0.2, turned out to be in better agreement with 
the general trend of the (stellar) metallicity-luminosity relation 
for dwarf irregular galaxies, indicating a "normal", although 
metal-poor, dlrr galaxy. 

The cold and warm ISM of SagDIG and its kinematics have 
been the subject of several investigations. Skillman et al. ( 1989 1 
obtained optical spectrophotometry of the most luminous Hn 
regions, and estimated a [O/H] abundance ~ 3% of the solar 
value, a result independently confirmed by the new measure- 
ments of Saviane et al. (I2002> . 

High-resolution, high sensitivity VLA observations of the 
SagDIG H i content have been obtained by Young & Lo (119971 
(see also Lo et al. [T993l . About 1 .2 x 10 7 M of (H + He) have 
been estimated by M02 using (m - M) = 25.14 and the Hi 
column density given by Young & Lo. SagDIG thus appears 
to have a relatively high mass fraction in the form of neutral 
gas, Mhi/Lb = 1.6, a value typical of the Hi content in dwarf 
irregular galaxies (see Mateo 1998J. The H i is distributed in an 



almost symmetric ring which was interpreted as produced by 
the combined effects of stellar winds and supernovae. 

This paper presents our new HST/ACS observations, 
clearly overcoming the limitations of previous ground-based 
observations and probing the young, intermediate age, and old 
stellar populations in SagDIG. The results of those observa- 
tions are presented in this paper. The paper is organized as fol- 
lows. In Sect. |3 we introduce our HST/ACS observations and 
data reduction. Section[3]presents the ACS color-magnitude di- 
agrams and their main features. In Sect.|4]we re-derive the dis- 
tance to SagDIG, the mean metallicity of its red giants stars, 
and, most importantly, provide evidence for a 10 Gyr old hori- 
zontal branch population. The recent star formation is then an- 
alyzed in Sect. |5] Section|6]compares the spatial distribution of 
the stars and the H i and discusses the relationship between the 
composite stellar populations of SagDIG and the morphology 
of its ISM. The principal results of our ACS study are finally 
summarized in Sect.0 




Fig. 1. A 3.1 x 2.4 arcmin true-color image of SagDIG derived 
from our ground-based ESO/NTT BVI observations (1" reso- 
lution). North is to the top, East to the left. The low resolution 
of NTT data highlights large-scale structures such as tails of 
massive hot stars and regions almost devoid of young stars. 



2. Observations and data reduction 

Images of SagDIG were taken using the ACS on board the 
HST, in the course of one visit (on August 18, 2003; program 
GO 9820) for a total of 3 orbits. The observed field covers 
the main body of the galaxy and is centered at I = 2L06, 
b = -16°28. Through the ACS-Wide Field Channel (WFC), 
we obtained 5x396 s images in the F606W filter, 5x419 s in 
F475W and 5x421 s in F814W. The series of 5 images in 
F475W and F606W were dithered so as to allow filling the 
gap between the CCDs, removal of cosmic rays, and optimiza- 
tion of the point spread function (PSF) and spatial sampling 
(Anderson & King 120001 1. No attention has been paid in the 
recovering the gap in F814W images (sub-arcsec dithering pat- 
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Fig. 2. True-color image of SagDIG derived from our ACS observations (see also http://heritage.stsci.edu/2004/31/index.htmli. 
A central, asymmetrical region of recent star formation, with faint extensions, is superimposed on a sheet of red giant stars 
belonging to the older stellar population of SagDIG. Superimposed are the Hi column density contours from Young & Lo 
(1997). Also plotted are the locations of the H 1 hole and H 11 region. 



tern), so as to avoid residuals in the geometric distortion cor- 
rection when second epoch F814W images (to be obtained in 
cycle 14) are used to derive the proper motion of observed stars. 

The instrument distortion was corrected using ACS/WFC 
geometrical distortion correction maps made by Anderson 
(2004, priv. comm.). All images in the same filter were then 
stacked to build a median image, with pixel oversampling by 
a factor of 2. These steps allowed us to improve the signal-to- 
noise ratio and the PSF sampling. Using these images, standard 
D AOPHOT/ ALL S TAR (Stetson [T9"g7t routines were employed 
to derive both the PSF magnitudes and their aperture correc- 



tions. A quadratically variable PSF was adopted to account 
for the PSF variations across the field of ACS (Krist I20U3I . 
The photometric calibration of the ACS/WFC filters follows 
Bedin et al. (j2004i. This zero-point photometric calibration is 
based on the Vega spectrum and follows the recipes given by 
Holtzmann et al. ( 1995 1 for WFPC2, using the most updated 
in-flight characterization of the camera (Sirianni et al. l2004l >. 

Photometric errors and completeness were estimated from 
extensive artificial stars experiments following procedures sim- 
ilar to those of M02. More than 50,000 stars were added to the 
stacked images in a range of magnitudes and colors match- 
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ing the distribution of SagDIG stars in the observed color- 
magnitude diagrams. The artificial stars were placed in a suit- 
able pattern of triangular meshes to avoid self-crowding (see 
Saviane et al. l2000t . The images with the artificial stars added 
were processed exactly in the same way as the original im- 
ages (in particular photometry and calibration were done using 
the same parameters) and the photometry catalog was matched 
with the input artificial star list. It should be noticed that this 
procedure is approximated since it does not fully reproduce 
the spatial variations of the PSF. Still, it should account for 
the typical errors and completeness of the data. The results for 
photometric errors and completeness are presented in the form 
of completeness contours in the color-magnitude diagrams (see 
below) and in Tabled 

Figure ^ shows a color-image of our ground-based study 
(M02), highlighting some large-scale patterns like tails and 
holes in the distribution of the young blue stars. Figure[2]shows 
a true color image of SagDIG derived from our HST/ACS 
data. The spatial resolution of ACS allows us to resolve both 
the young and old component in SagDIG. It also confirms the 
large-scale patterns seen in ground-based images. 

3. Color-magnitude diagrams 

The ACS color-magnitude diagram of SagDIG reveals for the 
first time its history of star formation (Fig. [3j. The high- 
resolution of HST/ACS allows overcoming the crowding diffi- 
culties, producing highly detailed CMDs where many distinc- 
tive features of the various stellar populations are seen for the 
first time. Figure[3]presents both the hzf606w> ( ot f475w -twf606w) 
and the mpsww, ('«F475w - '«F8i4w) diagrams, which comple- 



Table 1. Photometric errors and completeness as derived from 
artificial star experiments. 



mag 


C F475W 


Cf475W 


f f606W 


Cf 606 W 


O" F8141V 


CfS14W 


20.00 


0.010 


1.000 


0.010 


1.000 


0.010 


1.000 


20.40 


0.010 


1.000 


0.010 


1.000 


0.010 


1.000 


20.80 


0.010 


1.000 


0.010 


1.000 


0.009 


0.999 


21.20 


0.010 


1.000 


0.010 


1.000 


0.009 


0.999 


21.60 


0.010 


1.000 


0.010 


1.000 


0.011 


0.999 


22.00 


0.010 


1.000 


0.010 


1.000 


0.013 


0.997 


22.40 


0.010 


1.000 


0.008 


1.000 


0.013 


0.996 


22.80 


0.010 


1.000 


0.012 


0.999 


0.016 


0.996 


23.20 


0.009 


1.000 


0.014 


1.000 


0.021 


0.994 


23.60 


0.012 


1.000 


0.016 


0.999 


0.025 


0.991 


24.00 


0.019 


1.000 


0.019 


0.997 


0.033 


0.988 


24.40 


0.027 


0.996 


0.023 


0.996 


0.043 


0.989 


24.80 


0.036 


0.993 


0.028 


0.994 


0.052 


0.990 


25.20 


0.046 


0.989 


0.036 


0.992 


0.066 


0.963 


25.60 


0.060 


0.988 


0.047 


0.991 


0.086 


0.894 


26.00 


0.077 


0.974 


0.059 


0.988 


0.114 


0.801 


26.40 


0.097 


0.932 


0.076 


0.978 


0.148 


0.715 


26.80 


0.125 


0.807 


0.102 


0.966 


0.186 


0.603 


27.20 


0.167 


0.677 


0.137 


0.924 


0.225 


0.384 


27.60 


0.211 


0.428 


0.175 


0.552 


0.263 


0.168 


28.00 


0.256 


0.182 


0.212 


0.141 


0.302 


0.047 


28.40 


0.301 


0.043 


0.249 


0.008 


0.341 


0.008 


28.80 


0.345 


0.010 


0.286 


0.009 


0.379 
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ment each other in highlighting different evolutionary phases. 
The CMDs presented in this section are limited to stars with 
absolute SHARP values below 0.3. 

Some foreground contamination is present in all diagrams. 
This is dramatic in the case of ground-based observations, yet 
it does not significantly affect the ACS data. Indeed, because of 
the distance of SagDIG, most of the galaxy stars are fainter than 
the main-sequence turnoff(s) of the thick disk stars, which are 
the main contributor to the foreground counts. In this respect, 
ACS really allows us to look through the Milky Way. 

As an aid to the interpretation of our CMDs, we 
have simulated the color-magnitude diagram of the Galactic 
stars along the line of sight of SagDIG using the 
Besancon Galaxy Model (Robin et al. 120031 see also 
http://bison.obs-besancon.fr/modele/!. This is shown in the 
upper right panel of Fig. [3] The foreground stars belong to 3 
main populations: (i) the Galactic disk, (ii) the thick disk, and 
(iii) the halo. Although the simulated diagram cannot be em- 
ployed to quantitatively reproduce the observed Galactic con- 
tribution, because of the different pass-bands, the overall simi- 
larity of the simulated and observed diagrams is helpful to in- 
terpret the observed CMD features. The vertical sequence at 
(otf475w - '«F606w) = 1.4 is clearly produced by the superpo- 
sition of the lower main sequences of Milky Way stars pro- 
jected throughout the line of sight. The secondary peak located 
at (B - V)- 0.5 corresponds to the overlapping of G field main- 
sequence turnoff stars along the line of sight. For example, stars 
at »ip606w — 20.0 + 0.2 and »?f475w ~ »iF606w— 0.55 would cor- 
respond to the turnoff area of Galactic stars located at ~ 8 kpc 
from us. These stars are quite important on their own, as they 
provide clues on the Galactic structure. A thorough analysis of 
the projected distribution of Milky Way stars will be possible 
at the completion of our program, when second-epoch data will 
enable us to analyze the proper-motion of all the observed pop- 
ulations. 

Turning again our attention to SagDIG, we note that the 
red giant branch (RGB) is now well sampled from its tip at 
wf814w~ 21.3 down to mp8i4w~ 26, i.e. well below the red 
clump of core He-burning stars. This RGB is accompanied by a 
parallel and slightly bluer sequence of red super-giants (RSG), 
whose most luminous stars are confused with the Galactic fore- 
ground. The red super-giants are better seen in Fig. |4] (left 
panel), where we plot the CMD of stars in the inner region of 
SagDIG. This diagram is dominated by the young stellar pop- 
ulation. On the other hand, the CMD of the outer regions (see 
right panel of Fig.0} shows, as expected, a decreasing fraction 
of young stars. 

One of the most noteworthy features of the ACS color- 
magnitude diagrams in Fig. [3] is the perfect separation of the 
H-burning main sequence and the sequence of blue He-burning 
stars (BHeB) in the blue-loops evolutionary phase. This is best 
seen in the OTf814W> ('«F475w - '«F8i4w) diagram of Fig. [3] Near 
'«F8i4w= 23, the color difference between the two populations 
reaches (f«F475w-'«F8i4w)= 0.35 mag. Besides the photometric 
precision of ACS, this result is obviously due to the wide color 
baseline chosen for our observations. At otfsmw- 25, the MS 
shows a clear bifurcation. A small component departs from the 
main sequence and proceeds obliquely, crossing the blue loops 
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Fig. 3. HST/ACS color-magnitude diagrams of SagDIG. The upper-left panel shows the mF606W, (otf475w - '«F606w) diagram 
for stars with \S HARP] < 0.3. The thick lines represent the 80% and 50% completeness levels as derived from artificial stars 
experiments. Mean error bars from the same experiments are plotted as a function of »Jf606w- The open symbols represent the 
carbon stars identified by Demers & Battinelli (2002, open circles) and Cook J 19871 open squares). The thin lines delimit the 
approximate location of the instability strip on the CMD. The upper-right panel shows a synthetic V, (B - V) diagram in the 
direction of SagDIG, calculated using the WEB interface of the Besangon Galactic model (Robin et al. 2003 ). Different symbols 
refer to the contributions of the Galactic thin disk (open circles), thick disk (small squares) and halo (crosses). The lower panel 
shows the otf814w, ('«F475w - '«F8i4w) diagram. 



at mF8i4w~ 24.0. This feature is indicative of a ~ 200 Myr old 
burst of star formation (see Sect. A similar feature is seen 
(although at a slightly lower resolution because of the narrower 
color baseline) in Sextans A (Dolphin et al. l2003i . 

We also notice a gap in luminosity between mp8i4w= 20.6 
and 21.1, in the distribution of blue super-giants centered at 



(otf475w - '«F8i4w- 0.5). Although this hint is based on a sta- 
tistically small sample of stars, we note that any quiescent star- 
forming phase in the past 100 Myr would not be detectable 
along the MS, since successively younger stars would overlap, 
and therefore fill any empty region along the main sequence. On 
the contrary, evolved blue or red super-giants do not overlap in 
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m F475W 



Fig. 4. Enlargements of the color-magnitude diagram illustrating the stellar population gradient in SagDIG. Left: The CMD of 
stars within the star forming regions, near the center of SagDIG. The curved line delimits the selection of MS stars used in 
calculating the star formation rate (see Sect. 5.2). Note the presence of the extended blue-loop feature and its red counterpart 
parallel to the RGB, indicating a significant population younger than 1 Gyr. This young population is less prominent in the 
right panel showing the h?f606w, 0«f475w - »?F606w) diagram for stars outside the star forming regions. The richness of the red 
clump and its complex structure are indicative of a predominant intermediate-age population formed over an extended epoch. 
The noticeable, although minor, extension to lower luminosity is reminiscent of the red horizontal branches seen in most dwarf 
spheroidal galaxies. The solid lines show the approximate location of the instability strip. 



luminosity, therefore preserving any signature of a quiescent 
phase. Thus, the blue super-giant gap between »?F8i4w= 20.6 
and 21.1 (at m^isw — '«F8i4w- 0.5) might trace a short quies- 
cent phase that we estimate (via isochrone fitting) took place 
some time between 30-60 Myr ago. 

Another important feature in the CMD of SagDIG is the 
presence of an extended asymptotic giant branch (AGB) which 
traces a conspicuous intermediate age population. The AGB in 
SagDIG is partly masked by the foreground vertical sequence 
centered at (»?F475w - /«F606w) =1-4, and becomes really evi- 
dent only when ACS photometric catalogs are cross-identified 
with the results of carbon star surveys. Figure [3] provides a 
cross-identification of our ACS data with the carbon star sam- 
ples of Cook (1983 open squares) and Demers & Battinelli 
(2002. open circles). This comparison confirms the composite 
nature of the SagDIG carbon stars sample. As suggested by 
Cook J1987L the C stars in SagDIG fall in two classes: (1) a 
luminous, intermediate-age population reaching ~ 2 mag in / 
above the RGB tip; and (2) a bluer, low-luminosity component, 
overlapping the upper 1 mag of the RGB, probably belonging 
to an older population. We mention a third possibility for the 
origin of the carbon stars identified in the SagDIG field. Of 
the 3 stars that coincide with the vertical foreground sequence 
(upper-left panel in Fig.[3j, the brightest one can be a luminous 



object on the extended AGB. However, the two fainter stars, 
quite separated from all other carbon stars, might belong to 
the Galactic foreground. These are possibly related to the pop- 
ulation of faint high-latitude Carbon stars found in the Sloan 
Digital Sky Survey (Downes et al. l2004t . The possibility that 
these two stars are Galactic low-luminosity dwarf Carbon stars 
will be better disentangled once our proper-motion analysis of 
SagDIG is completed. 

The color-magnitude diagrams in Fig.[3]also show a promi- 
nent red clump (RC) of core He burning stars, with a complex, 
elongated shape. This structure (which will be quantitatively 
analyzed along with the star formation history of the galaxy 
in a separate paper - Held et al. 2005, in prep.) indicates that 
the bulk of SagDIG stars formed at intermediate ages (^ 1 
Gyr), in a relatively extended star formation episode. The ex- 
tension of the RC at higher luminosities, sometimes known as 
the "vertical clump" (VC), marks the evolution of younger stars 
(~ 0.35 - 0.63 Gyr) in the "blue loop" core He-burning phase. 

In addition to these components, we note a fainter exten- 
sion at mp606w= 25.70 + 0.05 (rapg^w ~ 25.2) identified with a 
red horizontal branch (HB). Again, spatial selection helped us 
to pinpoint specific stellar populations. The presence of an old 
red HB is better seen in Fig.0](right panel), where stars located 
outside the regions of active star formation have been selected. 
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Detecting the presence of a truly old population in SagDIG was 
one of the primary goals of this study and, indeed, this red HB 
implies that SagDIG hosts an old f ^ 10 Gyr) stellar population. 
Together, Figs.[3]and gjmay also suggest that a search for vari- 
able stars, in particular the old (> 10 Gyr) RRLyrae variables, 
might be worthwhile, even though the number of RRLyrae 
stars we expect from Fig.|4]is not high. In these figures, the blue 
and red edges of the instability strip converted from the theoret- 
ical HR diagram (Chiosi et al. lTffHl to the ACS/WFC system, 
have been over-plotted on the CMD. The "gaps" noticed along 
the HB [(mp475w - '«F606w)= 0.3-0.4, mp606w- 25.8] as well 
as along the "vertical clump" [(mp475w - '«F606w)= 0.35-0.45, 
'«F606\v- 24.8] show an almost perfect coincidence with the in- 
stability strip. In particular, the red edge appears to agree very 
well with a truncation in the red HB and blue-loop sequence. 
Indeed, short-period variable stars have been found in other 
dwarf irregular galaxies (e.g., Dolphin et al.f2002;i Skillman et 
al. 120031 Baldacci et al. 120041 . Unfortunately, a similar search 
cannot be conducted using our present observations because of 
their limited 3 hours baseline. Selection of candidate variable 
stars is likely to become possible once second-epoch observa- 
tions are available. 

4. Reddening, metallicity, and distance 

4.1. Reddening 

The range of E(B - V) values for SagDIG in the literature 
is not particularly large, spanning the interval between 0.06 
(Lee & Kim l2000l up to E(B - V) = 0.19, derived in Saviane 
et al. d2002t from the Balmer decrement in an Hn region. 
Nonetheless, this range may significantly affect our conclu- 
sions on the stellar metallicity in SagDIG. Momany et al. 
( 2002} derived E(B - V) = 0.07 from the colors of foreground 
stars at the MS turnoff, while the Schlegel et al. d!998> extinc- 
tion maps indicate a foreground reddening E(B - V) — 0.12. 
These two values are likely to bracket the foreground redden- 
ing toward SagDIG. 

In addition, some degree of internal reddening has been 
suggested by Momany et al. d2002t to be associated with the 
star-forming regions. Accordingly, a single value of redden- 
ing for all stellar populations in SagDIG may not be appropri- 
ate. A comparison of the ACS color-magnitude diagrams with 
de-reddened isochrones appears to confirm this earlier sugges- 
tion (Fig|51 upper panel). In order to test the reddening value 
needed to fit the young main sequence, a set of 63 Myr old 
isochrones with different metallicities are over-plotted on the 
SagDIG color-magnitude diagram in the upper panel of Fig. [5] 
The isochrones provided by L. Girardi are from the Padua stel- 
lar evolution models, calculated in the same way as Girardi et 
al. d2002t . using the ACS transformation curves by Sirianni et 
al. d2004t . These isochrones assume Y= 0.230. The conver- 
sion from E(B - V) to the extinction in the ACS/WFC bands 
was taken from Bedin et al. d2004> . 

In the upper panel of Fig.[5]we over-plot 4 isochrones with 
Z = 0.0001, 0.0004, 0.001 and 0.004 (respectively from left to 
right) shifted by E(B - V) = 0.12. One sees that assuming this 
reddening value even the un-plausible hypothesis of metal-rich, 
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Fig. 5. The upper panel shows theoretical isochrones in the 
ACS-Vega system on the m F gi4w, («f475w - twfshw) di- 
agram. Solid isochrones have been shifted using a (m - 
M) a = 25.14 and E(B - V) = 0.12, whereas the long-dashed 
isochrone has been shifted using E(B - V) = 0.20. Solid 
isochrones, and from left to right, are with Z = 0.0001, 
0.0004, 0.001 and 0.004 (corresponding to [Fe/H]=-2.36, 
-1.74, -1.33 and -0.72). The long-dashed isochrone is with 
Z = 0.0004. All isochrones refer to an age of 63 Myr. 
The lower panel shows the mp8i4Wj ('«F606w _ '«F8i4w) color- 
magnitude diagram of SagDIG with superimposed the fidu- 
cial lines of the Galactic globular clusters NGC 6752 (right) 
and NGC 6341 (left), for which we used the reddening and 
distance moduli given in the Harris d!996> on-line catalog: 
http://physun.physics.mcmaster.ca/~harris/mwgc.dat as up- 
dated in February 2003. 



47 Tuc-like, isochrone is still bluer than the bulk of the SagDIG 
main sequence. On the contrary, a reasonable Z = 0.0004 (as 
derived from the [O/H] of the brightest Hn region, Saviane et 
al. l2002t isochrone shifted by E(B — V) = 0.20 shows an excel- 
lent fit to the SagDIG main sequence. 

A similar comparison with isochrones cannot be employed 
for the RGB colors because of the well-known uncertainties in 
the color-metallicity calibration of theoretical isochrones for 
red giant stars (e.g., Westera et al. I2002i . However, we can 
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Fig. 6. The 2-cr selected sample of RGB stars in SagDIG, plot- 
ted on the Mp8i4w vs. (mp606w - '«F8i4w)o diagram. On the left 
panel, the data are shown along with RGB fiducial lines of the 
Galactic globular clusters NGC 6341, NGC 6752, NGC 104. A 
grid of analytic RGB templates are superposed to the data in 
the right panel (see text for details). 

compare the data for SagDIG with the RGB fiducial lines of 
Galactic globular clusters to check the foreground reddening 
of RGB stars. The lower panel of Fig. [5] over-plots the RGB 
fiducial lines of two Galactic globular clusters, NGC 6341 and 
NGC 6752, onto the color-magnitude diagram of SagDIG. The 
clusters were observed with the same instrumental setup (GO 
9453, RI. Brown; see Brown et al. I2003l > and re-calibrated 
by Bedin et al. (I2004> with the same methods as used here 
for SagDIG. Both NGC 6341 ([Fe/H] = -2.3) and NGC 6752 
([Fe/H] = -1.54) have been corrected to the distance of 
SagDIG using the reddening and distance moduli given in the 
Harris J1996i catalog, and shifted using E(B - V) = 0.12 and 
(m - M) = 25.14 to match the SagDIG diagram. While this 
comparison does not provide itself an independent measure- 
ment of reddening, it gives no indication of any internal red- 
dening outside of star-forming regions. Therefore in the case 
of RGB stars we adopt the Schlegel's et al. J1998i value as our 
estimate of the total reddening. In light of the results obtained 
in M02, we regard this reddening value as an upper limit and 
associate to it an uncertainty of +0.05 that brackets the low 
reddening value adopted in M02. 

4.2. Metallicity of the red giant stars 

The mean stellar metallicity of SagDIG was estimated by an- 
alyzing a sample of red giant stars selected in regions out- 
side the star-forming complexes. This was done using the 
Saviane et al. d2000> method, extended here to the F606W 



and F814W ACS pass-bands. The method is based on the 
construction of a family of hyperbolae in the plane M F gi4 W , 
( m F606w _ m F8i4w)o suitable to fit the observed RGB fidu- 
cial lines of template globular clusters in a wide range of 
metallicities. The present implementation makes use of the 
empirical RGB fiducial lines of 3 Galactic globular cluster 
(NGC 104, NGC 6752, and NGC 6341) observed with ACS. 
This set of globular clusters was originally observed by Brown 
et al. J2003t to constrain the metallicity of the M31 halo, and 
reanalyzed by Bedin et al. (|2004i. The adopted metallicity of 
the 3 clusters is [Fe/H] =-0.76, -1.56 and -2.28, respectively, 
as given in the Harris J1996l > catalog. The fiducial lines have 
been transformed into the Mp8i4w, ('«f606w - '«F8i4w)o plane 
by adopting (1) (m - M) v = 13.37, 13.13, and 14.64; and (2) 
E(B - V) = 0.04, 0.04, and 0.02 respectively. 

The application of this method to SagDIG in the Mpsuw, 
(?«F606w - f«F8i4w)o plane is illustrated in Fig.|6] where the ab- 
solute magnitude and reddening-corrected colors of individual 
stars in SagDIG are compared with the RGB ridge lines of the 
3 Galactic clusters (left panel), and the grid of analytic RGB 
lines (right panel). The observed magnitude and colors have 
been transformed into the Mpsnw, ('«F606w _ '«F8i4w)o plane by 
assuming E(B - V) = 0.12 and (m - M) = 25.14 (see below). 
The deepth of our ACS imaging allows a much better identifi- 
cation of the RGB with respect to ground-based observations. 
To prevent inclusion of young red super-giants contaminating 
the blue side of the RGB, we applied a 2cr rejection in color 
about a fiducial line drawn by eye. 

An interpolation over color provides a metallicity estimate 
for each star. The [Fe/H] values of the selected RGB sample 
span from [Fe/H] = -2.6 to -0.8. The resulting mean metallic- 
ity (for stars brighter than M F8 i4 W = -3.0) is <[Fe/H]> = -2.21 
with an r.m.s. of o~ [Fe/H] = 0.18. By repeating the calcula- 
tion with a lower reddening, E(B - V) = 0.07 (as in M02) 
we would obtain ([Fe/H]) = -1.89 + 0.17. The value mea- 
sured by Momany et al. d2002> was <[Fe/H]> = -2.10. Given 
the uncertainties and the different pass-bands and calibrations, 
we regard these results as perfectly consistent. According to 
the new [Fe/H] measurement, SagDIG is confirmed to be very 
metal-poor yet "normal" among dwarf irregular galaxies. We 
note that this metallicity estimate does not take into account 
the fact that the RGB stars in SagDIG are younger, on average, 
than the red giant stars in globular clusters. However, the ef- 
fect is quite small for metal-poor stars (see Rizzi et al. l2003> . 
We will implement this correction after discussing the SFH of 
SagDIG (E. V. Held et al. 2005, in prep.). 

4.3. Distance 

The excellent match between the RGB tip of SagDIG and that 
of GGCs obtained by assuming the distance derived by M02 
suggests that the new ACS photometry is consistent with our 
previous distance estimate. Given the superior RGB obtained 
from ACS data, it appears worthwhile to verify the distance 
to SagDIG based on the magnitude of the RGB tip. Figure 
shows the luminosity function of RGB stars in SagDIG to- 
gether with our measurement of the tip. The RGB selection 
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Fig. 7. The upper panel shows a comparison between the ACS 
'«F8i4w data and /ntt photometry from Momany et al. J2002> . 
Stars plotted with heavy symbols are a 3<x-clipped subsample 
yielding a mean shift of -0.03 + 0.03 mag. The lower panel 
shows the luminosity function of SagDIG RGB sample se- 
lected outside the star-forming regions. The thin line repre- 
sents the best fit with an error-convolved step function used 
to measure the RGB cutoff. An arrow marks the RGB tip at 
m F814w (RGBT) = 21.23 ± 0.05. 



is the same as used for measuring the metallicity. The RGB tip 
was found at »ifsi4w(RGBT) — 21 .23 + 0.05 by fitting an error- 
convolved step function as in M02. Since the color term for 
the filter F814W is small in the color range of RGB stars (see 
Sirianni et al. 2004). we have been able to compare the new 
F814W magnitudes with the ground-based /-band photometry 
of M02. The result is shown in Fig.0 upper panel. The median 
difference is very small, /ntt - '«F8i4w = +0.04 + 0.03. By ap- 
plying this correction, we obtain / t ; p (ACS) = 21 .27 + 0.05. This 
is only 0.03 mag brighter than the value measured in Momany 
etal. f2W2\ . hr.RT = 21.30. 

Assuming the new tip determination at / = 21.27 + 0.05 
and a reddening of E{B - V) = 0.12 + 0.05 (i.e. A/ = 0.21 
comprising almost all reddening values reported in the litera- 
ture) and adopting M^ GBT = -4.04 + 0.05 for metal-poor sys- 
tems like SagDIG (see Da Costa & Armandroff l 19901 Lee et 
al. 119931 Bellazzini et al. 12001b . the resulting distance modulus 
as derived from the calibration of the RGB tip is (m - M) a - 
25.10+0.1 1. The distance error includes all uncertainties on the 
RGB tip determination and absolute calibration and reddening 
through standard error propagation. This basically confirms our 
previousM02 distance estimate, (m-M) = 25.14 + 0.18, - 1 .1 
Mpc. The 0.04 smaller distance modulus obtained using the 
ACS data is to be attributed to a slightly brighter (0.03 mag) 
RGB tip determination, the assumed M^ GBT , and 0.05 higher 
E(B - V). We note however, that calibrations of the absolute / 
magnitude of the RGB tip are done in standard ground-based 



/ Cousins filter. Since the transformation from the ACS mag- 
nitudes to standard / is quite uncertain, we retain our ground- 
based distance determination of SagDIG, and this is adopted 
for the rest of the paper. 

4.4. The Horizontal Branch 

Using this distance, Fig. [8] provides convincing evidence that 
the fainter extension of the RC identified in SagDIG is indeed a 
red horizontal branch. In order to estimate the magnitude level 
of the horizontal branch, we extract from the CMD stars in the 
color range (mp475w - '«F606w) = 0.40 - 0.55. The magnitude 
distribution of these (mostly core He-burning) stars is shown in 
Fig.[8](upper-right panel) as a histogram with bin size 0.5 mag. 
A multi-bin luminosity function is plotted to reduce the effects 
of the choice of the bins. The red horizontal branch is clearly 
detected at f7?F606w = 25.70 + 0.05. This is confirmed when 
we over-plot, in the upper-left panel, a Zero- Age HB (ZAHB) 
model of [Fe/H]^ -2.1 from Pietrinferni et al. (120041 . having 
assumed E(B - V) = 0.12 and (m - M) Q = 25.14, that perfectly 
fit the estimated level of HB in SagDIG. 

Further evidence for the correct identification of the HB is 
given in the lower panels of Fig. [H] where we compare the HB 
morphology in SagDIG with (i) the HB of an old, metal-poor 
Galactic globular cluster, namely NGC2419, and (ii) the ACS 
observations of variable stars in M31 (Brown et al. l2004> . 

NGC2419 is one of the few clusters matching the SagDIG 
metallicity and observed in the same ACS filters (GO 9666, RI. 
Gilliland). Also plotted is a sketch of the instability strip. For 
the HB stars in NGC 2419 bluer than »tF475w -'«f814w= 0.5 we 
only plot a fiducial line because of the overlap with the young 
MS stars in SagDIG. This comparison shows the consistency 
of the location the HB level in SagDIG with that observed in 
the oldest globular clusters. On the other hand, it lends support 
to the adopted distance modulus and reddening. 

In addition, in Fig.|S]our data are compared with ACS ob- 
servations of RR Lyrae stars in M31 (Brown et al. [2004). In 
this case, we have arbitrarily adjusted the M31 data to match 
the red-HB of SagDIG. This comparison suggests where one 
should search for RR Lyrae variable stars. 

5. The recent star formation 

5.1. Blue loop stars: age and metallicity 

Theoretical isochrones can now be used to characterize the 
ages of the MS and blue HeB stars. Since the metallicity of 
SagDIG is known to be in the metal-poor regime, we em- 
ployed isochrones with Z = 0.0001 and 0.0004 ([Fe/H]=-2.31 
and - 1 .73, respectively) for our comparison with observational 
data. Very metal-poor isochrones (Z = 0.0001) with ages 
31.5, 100, 200, and 630 Myr are over-plotted on the mpsnw, 
(otf 4 75w-'Wf814w) CMD of SagDIG in the upper panel of Fig. [9] 
while isochrones with Z = 0.0004 and ages of 3 1 .5, 89, 200 and 
630 Myr are shown in the lower panel. The MS and blue loops 
are quite well reproduced by both sets of isochrones from the 
brightest stars (mp8i4w= 20.5) down to the "vertical clump". 
However, the excursion of the blue loops in the isochrones, 




Fig. 8. Enlargements of the CMDs illustrating the detection of an old horizontal branch in SagDIG. The upper-left panel highlights 
the color interval used for measuring the magnitude distribution of He-burning stars, along with a [Fe/H]^ -2.1 ZAHB model 
(assuming Y= 0.245) from Pietrinferni et al. (12004b . In the upper-right panel we plot the histogram and a multi-bin, smoothed 
luminosity function of stars selected in upper-left panel. A line at OTF606W ~ 26.70 marks the level of the old red-HB. In the 
lower-left panel, we have superimposed on the SagDIG data the fiducial loci for the globular cluster NGC 2419 (Momany et al. 
2005, in prep.). These data were fully calibrated in the same way as for the SagDIG stars. The HB stars in NGC 2419 in the color 
interval 0.4 < mp475w — '«F8i4w < 1.4 are shown as open squares. The Brown et al. (I2004> data for RR Lyrae variable stars in 
M31 are shown in the lower-Right panel {crosses: RRab; triangles: RRc). Also plotted is the fiducial line for static HB stars in 
M 31. The Brown's et al. data have been arbitrarily shifted to match the red-HB in SagDIG. 



which is strongly correlated to metallicity, favors a less extreme 
metal-poor composition (Z = 0.0004, or [Fe/H] ~ —1.7) for the 
100-600 Myr old stars in SagDIG. The details of stellar evolu- 
tion of HeB intermediate-mass stars depend on many uncertain 
physical parameters (such as the amount of convective over- 
shooting, stellar rotation, and mass loss), so evolutionary tracks 
cannot be safely used to measure the metallicity of young HeB 
stars. However, the Padua model isochrones are known to pro- 
vide a good fit to systems with known age. Therefore, assum- 
ing that metallicity is the main driver of the morphology of the 
blue loops (but having in mind all the systematic uncertainties 
on the model physics), this comparison would suggest a mod- 
est chemical enrichment (~ 0.4 dex) in stars since the main 



star formation episode from which most of the RGB stars orig- 
inated. 

5.2. Star formation rate from main sequence stars 

While a full analysis of the star formation history of SagDIG 
based on synthetic color-magnitude diagrams will be presented 
elsewhere, we give here a first account of the mean star forma- 
tion rate (SFR) in the last 0.6 Gyr. The SFR was computed as 
in Saviane et al. (J2004J hereafter SHR04), using the equation: 

log h = (a C sp - ffssp) AM V - log At ( 1 ) 
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Fig. 9. Model isochrones from Girardi et al. superimposed on 
the mp8i4w, ('«F475w - '«F8i4w) diagram of SagDIG. Upper 
panel: Z = 0.0001 isochrones with ages of 31.5, 100, 200, 630 
Myr. Lower panel: Z = 0.0004 isochrones with ages 31.5, 89, 
200, 630 Myr. 



where assp and acsp are the slopes of the luminosity function 
(LF) of a simple (isochronic) and a composite stellar popula- 
tion, respectively. AMy and At are the luminosity and age in- 
terval over which the mean SFR is calculated. All the parame- 
ters in the equation are measured on the V main sequence LF, 
except for At which is obtained from the theoretical isochrones. 

The main-sequence LF of SagDIG was obtained by se- 
lecting MS stars as indicated in Fig. |4] and then construct- 
ing the My histogram with a bin size of 0.5 mag (Fig. I10l > 
1 . The absolute magnitudes were computed as usual from the 
relation My = V- (m - M) a -R v E(B - V), adopting a cor- 
rected distance modulus (m - M) = 25.14 + 0.18, a redden- 
ing E(B - V) = 0.19 + 0.04 (appropriate for the MS stars, see 
Sect.|tj}, and a ratio A v / E(B-V)= R v = 3.1. The mean slope 
of the LF of SagDIG, a C sp = 0.486 + 0.022, was calculated by 

1 We converted the mF606w photometry of MS stars into V magni- 
tude by computing a mean difference (Vntt ~~ «?F606w) = +0.05 mag 
(Momany et al.|2002). This procedure is acceptable for MS stars since 
we are considering a narrow color range around (B - V) ^ 0, while in 
general the mismatch with the Johnson V band is too large to allow a 
meaningful transformation. For example, for the Galactic foreground 
sequence at (m V475W - m F60 6w) - 1-40, <Vntt - "iF606w> = -0.50. 
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Fig. 10. The Mf606w luminosity distribution of MS stars in 
SagDIG (shaded histogram). The thick solid line is our best- 
fitting model, which can be regarded as the superposition of 
luminosity functions of many SSPs (thin solid lines). The LF 
of the young LMC cluster NGC 2004 is shown as a dashed his- 
togram. 



fitting a straight line to the LF for V < 26. This ensures that 
stars on the sub-giant branch do not contaminate the MS sam- 
ple. Given this limit, the useful LF interval is AMy = 5 mag. 
In order to calculate the corresponding age interval, we used 
the interpolation of the MS termination against age given in 
SHR04. The luminosity range of the LF corresponds to age lim- 
its from 676.7 Myr ago to 20.4 Myr ago, implying Af = 656.2 
Myr. A slope qtssp = 0.32 was adopted for simple stellar popu- 
lations, as in SHR04. 

According to Eq. \T\ we find a SFR h - 1.032 x 
10~ 8 SSPyr 1 , normalized to a simple stellar population. This 
means that in the quoted time interval (~ 0.6 Gyr), 6.8 
SSPs were created. Finally, based on the zero-point differ- 
ence between the LF of our template (the young LMC clus- 
ter NGC 2004) and a single SSP, we find that one SSP here 
contains 1.32 clusters like NGC 2004, or 2.65 x 1O 4 M (see 
SHR04 for further details on the method). Therefore the mean 
star-formation rate in the quoted interval is 2.65 x 10 4 x h — 
2.73 x 10- 4 M o yr -1 . 

Karachentsev et al. JT9991 KAM99) computed the SFR of 
SagDIG through a simulation of the observed color-magnitude 
diagram, and found SFR=13 + 5 x lO^M^r 1 , 6.6 + 0.8 x 
lO~ 4 M yr _1 , and 1.3 + 0.1 x lO _4 M yr _1 in the age ranges 
- 0.05 Gyr, 0.05 - 0.2 Gyr, and 0.2 - 15 Gyr, respectively. In 
order to compare these values with our rate, we computed the 
time-weighted SFR of KAM99 within our time interval (0 - 
0.677 Gyr), which is 3.3 x lO^Moyr -1 , or ~ 20% larger than 
the SFR computed above. Now KAM99 point out that the age 
range 0.2 - 1 Gyr is poorly sampled by their CMD, so our value 
is perfectly compatible with their SFRs. A further explanation 
for the reported difference is that KAM99 computed their SFR 
over an area which is larger than our ACS field. 

If we restrict our SFR computation to the - 0.24 Gyr in- 
terval (i.e. cutting the LF at V — 24.5), we obtain a rate sim- 
ilar to the previous one (2.6 x 10 _4 M Q yr _1 ), while KAM99 
find 6.6 x 10~ 4 M o yr _1 . One possibility is that the larger SFR 
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of KAM99 is due to an incorrect evaluation of the MS/blue- 
loop ratio. This would not be surprising, given the large errors 
affecting ground-based photometry, that do not allow a clear 
separation between the two sequences (as it is the case in the 
ACS data). We also note that KAM99, using the Her flux mea- 
sured by Strobel et al. (1991), i.e. 28.1 x 1(T 15 ergcirrV 1 , 
and applying the method of Aparicio et al. (2000), find a cur- 
rent SFR=1.5 x 10 -4 Moyr~', much closer to that which we find 
here. In summary, our method predicts that the recent SFR re- 
mained almost constant during the last few hundreds of Myrs. 
The SFH of SagDIG will be further considered with CMD sim- 
ulations of the present ACS data in a future paper. 

6. The spatial distribution of gas and stars 

The interplay between stars and the surrounding interstellar 
matter (ISM) in the form of neutral, ionized, and molecular 
gas, is a key aspect in understanding what triggers star for- 
mation and drives the wide range of properties observed in 
dwarf galaxies. To investigate this connection, we analyze here 
the distributions of stars of different ages in SagDIG from our 
HST/ACS observations and compare them with the distribution 
of the ISM. 

6.1. The distribution of the gas 

We begin our discussion by summarizing here the properties of 
the ISM in SagDIG. The observations of Hi in dwarf galaxies 
often decompose the neutral H into two components (see, e.g., 
Young et al. U997l and references therein). In the particular case 
of SagDIG, a warm (broader and more turbulent) Hi compo- 
nent having a velocity dispersion cr — 10 km s _1 was observed, 
along with a cold (narrower and more quiescent) phase with 
cr - 5 km s _1 (Young & Lo ll997t . The broad component is dis- 
tributed all over the galaxy, whereas the narrower component is 
concentrated in few H i clumps of about 8 x 10 5 M Q (6 clumps 
are listed in Table 1 of Young & Lo 1997). A prominent Hi 
clump (with a column density 1 . 1 x 10 21 cm -2 ) is nearly coinci- 
dent with the biggest, and probably only, H n region in SagDIG. 
This H ii region is 100 pc across, and 100 pc away from the H i 
clump center (Strobel et al. I1991i . These Hi clumps form an 
almost complete asymmetric ring (or shell) structure surround- 
ing a region of lower H i column density, hereafter referred to 
as the "Hi hole". The roundness of the outer Hi distribution 
led Young & Lo to suggest an almost face-on orientation of 
SagDIG. The diameter of the ring is nearly 4 arcmin (1.3 kpc), 
while the largest extent of the H i at a column density of 5 x 10 18 
ctrT 2 is 3.3 x 3.0 kpc. 

Young & Lo J1997I found very little systematic rotation in 
SagDIG, i.e. any rotation velocity is of the order of, or smaller 
than, the intrinsic velocity dispersion. However, some degree 
of symmetry in the velocity field led the authors to suggest the 
presence of some large-scale organized motion. 

Young & Lo i 19971 1 also suggested that the cold dense 
atomic phase tells us where conditions are right for the forma- 
tion of molecular gas and stars. The sizes of the Hi cold com- 
ponents imply virial masses which are a few times larger than 
the H i mass alone; the difference may be caused by molecular 



gas. Thus, the cold H i component may trace the raw ingredient 
necessary for any star formation activity. Indeed, it is found in 
clumps of about the same size as the star formation regions. 
In dwarf irregulars, star formation typically occurs in regions 
of 100-200 pc in size, with each episode lasting for about 100 
Myr and then moving around (see, e.g., Dohm-Palmer et al. 
I1997ll998l for the case of Sextans A). 

6.2. Spatial distribution of young and old stellar 
populations 

The unique photometric deepness and precision of our ACS 
data allowed us not only to disentangle stars in different evo- 
lutionary phases (e.g., main-sequence stars against blue super- 
giants) but also to assign reliable age estimates to them from a 
comparison with the Padua isochrones on the WFC/ACS sys- 
tem (L. Girardi 2004, priv. comm.). Figure ITTI shows the dis- 
tribution of stars in different age intervals superimposed on the 
H i contour map of Young & Lo J1997I . kindly provided by 
the authors. Astrometric calibration of the radio map and our 
optical F606W reference image was performed independently 
and carefully checked. In all panels, the HST/ACS footprint is 
superposed on the H i map so as to mark the boundaries of our 
field of view. We note that, in order to prevent distribution bi- 
ases caused by the spatial variations in the PSF of WFC/ACS, 
no selection based on image shape parameters was applied to 
the photometry catalogs when plotting these spatial distribu- 
tions. For instance, we found that spurious distributions may 
originate from selection on the \S HARP] parameter due to the 
spatially varying image quality. 

The location of the youngest stars in SagDIG, i.e. the main- 
sequence stars with age less than 70 Myr (mpeoew^ 24.0), is 
plotted in Panel A of Fig.[TT] The MS and BSG stars were se- 
lected on the left and right side, respectively, of the separation 
line outlined in Fig. 0] The most evident feature is the star- 
forming region located in between the H i clumps ii and Hi on 
the eastern part of SagDIG. The angular size of this complex 
is ~ 30", or 155 pc, in diameter. Very recent star formation is 
also seen in the central part of the optical galaxy, on the rim 
of the H i "hole", in between the H i clumps iv and v. A very 
small complex, or a star cluster, with only a few very massive 
stars is seen close to the SE hydrogen clump (/). This clus- 
ter is of particular interest being associated with the only H n 
region in SagDIG. Panel B plots the distribution of MS stars 
having age between 70 and 160 Myr (24.0 <m F 606w^ 25.0). 
The stellar distribution starts to be more widespread, with a 
tail extending out to the Hi clump (v) on the western side. 
The blue-loop super-giants in the same age range 70-160 Myr 
('«F606w ^ 23.0, < m F 475w -'«F606w ^ 0.2), shown in Panel C, 
are consistent with this pattern, although their number is very 
small. Panel D of Fig. II II summarizes the distribution of the 
most recent star formation in SagDIG (age < 160 Myr), show- 
ing its nearly circular symmetry, with a prominent eastern star- 
forming complex and a peculiar tail of massive stars reaching 
the westernmost hydrogen cloud. 

The distribution of older MS stars (25.0 < mp606w ^ 26.5) 
and blue super-giants (22.5 < mF606w ^ 24.8, mp475w _ 
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Fig. 11. Spatial distribution of MS and BSG stars of different ages superimposed on the H i contours from Young & Lo d!997> . 
North is to the top and East to the left. Also plotted is the HST/ACS footprint. Panel A: MS stars with WF606W < 24.0 and age less 
than 70 Myr. Panel B: MS stars with 24.0 < mp606w ^ 25.0 and age between 70 and 160 Myr. Panel C: BSG stars with rap606w ^ 
23.0, age 70-160 Myr. Panel D: all MS and BSG stars younger than 160 Myr. Panel E: MS stars with 25.0 < m F6 06w < 26.5, age 
interval 160-700 Myr. Panel F: BSG stars 22.5 < m Fmw < 24.8, age between 160 and 700 Myr. 



7WF606W < 0.4; age between 160 and 700 Myr) in shown in 
Panels E and F. The eastern star-forming region, which is main 
site of very recent star formation, is also present for stars older 
than 160 Myr. In this age range, however, a complex on the 
western side of H i cloud iv becomes also prominent. The latter 
appears to have been the most active star-forming region more 
than a few hundred million years ago. 

Figure ^2 gives us a clear picture of the relationship be- 
tween the Hi and the recent star formation (< 700 Myr) in 
SagDIG. It confirms the earlier finding by Young & Lo d!997i 



that recent star formation sites are always near (although not 
coincident) with the H i column density peaks. This is best seen 
for the eastern and western SF complexes, as well as the small 
association coincident with the H n region. 

The spatial distribution of relatively old stars, including 
main-sequence stars older than 0.7 Gyr, He-burning stars on 
the red clump and the red-HB, and RGB stars, is shown in 
Fig. [21 Panel A shows the distribution of bright RC stars 
(24.8 < »?F606w < 25.4, 0.5 < mp475w - ?«F606w < 0.7). A com- 
parison with isochrones indicates for these stars an age interval 
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Fig. 12. Spatial distribution of relatively old stars. Panel A shows the distribution of bright RC stars with 24.8 < mF606w < 25.4 
and 0.5 < (»zf475w - '«F606w) < 0.7, corresponding to an age interval between 0.3 and 2.0 Gyr. Panel B shows the distribution of 
presumably old red-HB stars, having 25.6 < mp606w ^ 25.8 and 0.4 < (mp475w - '«F606w) < 0.6. Panel C shows the distribution 
of RGB stars brighter than mp606w = 24.5. Red giant stars are at least 1 Gyr old. For comparison, Panel D shows the distribution 
of MS stars with age between 0.7 and 2.0 Gyr (26.5 < m F60 6w ^ 27.5). 



between 300 Myr and 2.0 Gyr. The distribution of the old red- 
HB stars (25.6 < m F606w < 25.8, 0.4 < m F47 5 W - m F60 6w ^ 0.6) 
is shown in Panel B. These stars are most likely older than 
9-10 Gyr. A similar plot is done for RGB stars brighter than 
'«F606w = 24.50 in Panel C. 

For the stars younger than 1 Gyr, a correlation with the dis- 
tribution of the gas and the star-forming regions can still be 
noticed. The red clump stars, for instance, appear to follow the 
the rims of H i clumps on the eastern side. Yet, there is no par- 
ticular hint of a H i hole imprint. This trend is not present any- 
more in RGB stars (older than 1 Gyr) which show some central 
concentration. For the genuinely old (red-HB) population the 
statistics is really poor - all we can say is that their distribution 
is rather uniform and extended. 



It is interesting to compare these distributions with the lo- 
cation of MS stars older than 0.7 Gyr (26.5 < mpeoew < 27.50, 
Panel D). The fainter limit roughly corresponds to our 50% in- 
completeness level, or an age of about 2.0 Gyr. The Hi hole 
region is smoothly populated, and the star distribution is - al- 
ready in this age interval - fairly uniform with an overall central 
concentration, and extended. There are no detectable voids in 
the distribution of stars. However, the distribution is clearly in- 
homogeneous, still showing the presence of the north-western 
clump, and a generally higher surface density of stars on the 
eastern region. This implies that the sites of recent star for- 
mation discussed for the younger stars, in particular the north- 
western site, are quite old being already in place 1-2 Gyr ago. 
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6.3. The H i Hole 

If the large Hi hole in SagDIG is a structure that has formed 
over time from an initially more uniform distribution we can 
attempt to address the question of how it formed using con- 
straints from our deep ACS imaging. Before we continue how- 
ever we should stress that the assumption that the hole formed 
from a uniform distribution of atomic gas may be unfounded. 
This point is discussed at the end of this section. 

A popular scenario for the formation of H i holes is as the 
result of the energy input to the ISM from massive stars in the 
form of stellar winds and more importantly supernovae. Several 
galaxies have been observed with H i holes, and in some cases 
these holes are observed to be expanding, although this is gen- 
erally observed only for small holes (refer to HoII, M33, LMC, 
SMC, M31, Sextans A). Opinion is divided, even for a given 
galaxy, as to whether the supernova-driven wind hypothesis is 
viable or not [see for example Rhode et al. J1999b and Stewart 
etal. (HHP]- 

If the H i hole in SagDIG formed from the combined ener- 
getic input of supernova explosions we might expect to detect 
the light from those stars which formed along with those that 
drove the expansion but that have not exploded as supernovae. 
We would expect these stars to be found in a cluster or clus- 
ters that should still be somewhere within the H i hole (it is ex- 
tremely unlikely that such clusters would have dispersed com- 
pletely during the formation time of the hole). Because the stel- 
lar surface brightness within the Hi hole in SagDIG is clearly 
very low we will make assumptions in this estimate so as to 
minimize the stellar luminosity which should remain. 

From the Hi properties of the Hi hole in SagDIG we can 
estimate the energy required to create it from the expression 
given by Heiles J1979i : 

£ E = 5.3xl0 43 «i 12 (^) 312 ^ 4 erg (2) 

where d is the diameter of the hole in parsecs, no is the hy- 
drogen particle density in ctrT 3 and cr v is the velocity disper- 
sion of the atomic gas (in the case of SagDIG no expansion 
of the Hi hole is detected). This is a lower limit if there is 
expansion in the plane of the galaxy which is hidden by the 
face-on orientation of SagDIG. The detailed H i observations of 
Young et al. have allowed all of these to be measured; they find, 
d = 1.3 kpc, <x v =6-11 kms~' and «o - 0.22 cirT 3 for the 
higher velocity dispersion component. Taking <x v = 10 kins -1 
we find for the hole in the SagDIG £ E ~ 1 -5 x 10 53 erg. 

Assuming that a SN supplies 10 51 erg (McCray & Kafatos 
119871 of mechanical energy to the ISM this energy require- 
ment implies that it was created by ~ 150 SNe. Assuming a 
normal (Salpeter) IMF we therefore normalize the stellar mass 
distribution so that there are 150 stars more massive than 7 M , 
which gives an IMF, <£ m (ra) = 2880 m~ 235 . 

For this IMF the luminosity of all the stars less massive than 
7 M Q will depend on the age of the population. Although this 
will fall with time the reduction in luminosity will be modest 
over the lifetime of the hole. For the moment we assume the 
hole has the same age as the lifetime of the lowest mass pro- 
genitor that undergoes a supernova, ~ 5 x 10 7 yr (Oey & Clarke 
[T997l but see below. 



We estimate a lower limit to the luminosity of the remaining 
stars just after the last supernova occurred by assuming that 
they are all still on the main sequence. Using a standard main 
sequence luminosity function and integrating from 0. 1 to 7 M 
we find a luminosity of 4500 L Q . In order to reduce this number 
we must invoke aging of the population, but stellar evolution 
models show that between 5 x 10 7 and 10 8 yr this amounts to 
less than a factor of 2. 

At the distance of SagDIG 4500 L corresponds to m v =s 
21. Within the Hi hole we see no evidence of any clusters 
brighter than ra v =* 26.0 and so if the star formation episode that 
formed the hole was confined to one or even several clusters 
these clusters must have completely dispersed on a timescale 
of 5 x 10 7 years. Perhaps the star formation that formed the 
hole was not confined to clusters but was more widespread? 
Although the total stellar light integrated within the Hi hole 
does exceed 21 mag this light is distributed smoothly with an 
intensity that increases gradually towards the optical center of 
the galaxy. The surface brightness is always below that in other 
regions of the galaxy. For the supernova hypothesis to be feasi- 
ble then, we require that the stars with ages in the approximate 
range 50-100 Myr are concentrated within the Hi hole region 
of the galaxy. In Fig. we see that the only stars in this age 
range within the H i hole are located at its NW edge and that 
there are more stars in this age range beyond the confines of 
the hole. Indeed there is no evidence of any age stellar popula- 
tion that is preferentially distributed within the Hi hole. 

Possible ways to reconcile the data with the supernova- 
driven wind hypothesis. 

- The cluster or clusters in which the stars formed have been 
completely dispersed within ~ 5 x 10 7 yr. For a stellar clus- 
ter to be so short-lived one would expect a dynamic dis- 
ruption (such as the passage through a spiral arm or tidal 
interaction) to be the cause. Neither of these seem to be 
present in SagDIG and so one would expect clusters to be 
long lived. 

- The hole is very old. As the H i hole shows no evidence for 
expansion it could have formed some time in the past and 
just be a relic of an ancient star formation event. However, 
the velocity dispersion of the H i which is observed to be 
10 kins 1 places a limit on how long structure can survive. 
In the case of the hole in SagDIG, with a radius of 650 pc 
this is ~ 6x 10 7 yr. Therefore, adding the time of formation, 
the stars that formed coevally with those that expanded the 
hole should not be older than ~ 10 8 yr. Stellar evolution 
models show that for the luminosity of a cluster to fall by 
a factor of 100 (from mag 21 to our detection limit of 26) 
requires ~ 2 x 10 9 yr. We cannot therefore reconcile our 
non-detection of clusters with the presence of the H i hole. 

- A very top-heavy IMF. The analysis in first part of this 
paper gives no support for anything but a normal IMF in 
SagDIG. 

- The energy imparted to the ISM by a supernova is roughly 
2 orders of magnitude above the value we have used of 
10 51 erg, or the energy required to produce the hole is lower 
by a similar amount. Both of these would require that theo- 
ries of supernovae be substantially revised. 
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In the case of SagDIG the depth of the ACS data allow us 
to conclude that the supernova driven wind hypothesis is an ex- 
tremely unlikely explanation for the shell structure in SagDIG. 
Other authors have come to similar conclusions for the forma- 
tion of Hi holes (Crosthwaite et al. 120001 Deul & den Hartog 
[T990I Rhode et al. lT999l but that in SagDIG is surely one of the 
most secure cases of a failure of the supernova wind hypothe- 
sis. 

There are some possible alternatives to the supernova- 
driven wind hypothesis. 

- One possibility is that the H 1 hole region represents a re- 
gion where the ISM is dominated by molecular rather than 
atomic gas. Although the presence of molecular gas may 
not necessarily result in a detection of CO emission (CO 
is a poor tracer of molecular gas in low metallicity envi- 
ronments, Taylor et al. 1998 Rubio et al. I2004> one would 
expect to see some evidence of recent star formation asso- 
ciated with the molecular gas. The absence of young stars 
in this region therefore argues against this interpretation. 

- Another possibility is that the large gaseous ring is the re- 
sult of a dynamical resonance in a non-axisymmetric poten- 
tial. For such rings to form however the gas surface density 
must exceed the critical density given by the Toomre stabil- 
ity criterion. We cannot calculate this without knowledge of 
the rotation curve and so for SagDIG this cannot be tested. 
However, the gas surface density rarely exceeds the criti- 
cal density in dwarf irregular galaxies in which the rotation 
curve has been measured (e.g. ESO215-G009, Warren et al. 

mug. 

- A mechanism that may apply to SagDIG and also to dwarfs 
containing a single large H 1 "hole" such as the Local Group 
galaxies Sextans A (van Dyk et al. 119981 and LeoA (Young 
& Lo 119961 but probably not to dwarfs like Ho II which 
contain many smaller holes could simply be that the Hi 
distribution was never uniform in the first place. The accre- 
tion of an intra-group H 1 cloud by a low mass dark matter 
halo (in which some star formation had taken place at ear- 
lier epochs) would result in a far from uniform H 1 distri- 
bution. The non-zero angular momentum of the cloud with 
respect to the halo would naturally lead to the formation 
of a ring of gas. The ring would form at the radius where 
the gravitational force balanced the orbital angular momen- 
tum. For atomic gas that has been recently acquired from 
an 'external' source there is little reason to expect any cor- 
relation between the stars and the gas in the resulting sys- 
tem. Interaction between accreted gas and pre-existing ISM 
could have triggered the most recent star forming event 
(Fig. ^2p ane l A). The Hi structure around the Phoenix 
dwarf mapped by Young & Lo (I1997> in which 2 large 
clouds of H 1 may be present, neither of which is coincident 
with the stellar distribution, is consistent with this scenario. 
Rhode et al. dl999> discuss several alternative explana- 
tions for the formation of large H 1 holes in Ho II, including 
gamma-ray bursts, impacts of high velocity clouds (HVCs), 
large scale turbulence and ionization by the intergalactic ra- 
diation field. These possibilities are all equally applicable 
to SagDIG. Apart from noting that SagDIG is very isolated 



and that there is no evidence of nearby H 1 clouds we can 
do little to test these ideas. 



7. Summary 

We have presented a study of the stellar populations in SagDIG 
to unprecedented deepness using observations obtained with 
the Advanced Camera for Surveys aboard HST Our HST/ACS 
photometry is deep enough to go fainter than the main- 
sequence turnoff of halo, thick and thin disk stars in the Milky 
Way. Indeed, using ACS we can observe SagDIG looking 
through the Galaxy. Our deep color-magnitude diagrams reach 
down to mp606w ~ 27, mp8i4w ~ 27, which allowed us to ob- 
serve the stellar populations of SagDIG to a surprising degree 
of detail. The main results are the following: 

- As a result of the photometric precision of ACS and our 
choice of a wide color baseline, ffiF475w -ffl F8i4w> the young 
main-sequence is traced down to mp606w ~ 26, well sepa- 
rated from the sequence of He burning blue super-giants on 
the blue-loop phase. This allowed us to estimate the mean 
star-formation rate in the last 0.6 Gyr. 

- There is a well defined RGB, indicating that the domi- 
nant stellar population if older than ~ 1 Gyr. From the 
mean color of RGB stars we infer a mean stellar metallicity 
[Fe/H] = -2.21 (uncorrected for age effects). 

- Previously known carbon stars are easily identified on the 
CMD. Most of them belong to an extended AGB typical 
of an intermediate age population. However, a few C stars 
located with high precision below the tip of the RGB may 
belong to an old population. 

- The distribution of He burning stars in the CMD is indica- 
tive of an extended period star formation, spanning a time 
interval between as early as 10 Gyr ago and ~ 100 Myr ago. 
Most of the stars populate a very rich "red clump", which 
is the locus of stars older than ~ 1 Gyr. The stars younger 
than 1 Gyr are more luminous than the RC, evolving on 
the blue-loop phase to form the sequence of blue super- 
giants, and the associated (less prominent) sequence of red 
super-giants, contaminating the blue side of the RGB. The 
luminosity distribution of the blue super-giants will be used 
in a separate paper to reconstruct the history of star forma- 
tion up to 1 Gyr ago. The wide color range of blue-loop 
stars suggests that even the recently formed stars are sub- 
stantially metal-poor, among the most metal-poor in dwarf 
irregular galaxies (Z a 0.0004). 

- The presence of a red HB, noticed as a fainter extension 
of the RC - proves that an early epoch of star formation 
took place in this dwarf irregular galaxy. This result is of 
foremost importance since it sets a lower limit to the first 
epoch of star formation in yet another actively star-forming 
dwarf galaxy. In this regard, the resulting scenario is one in 
which all Local Group galaxies show to have had an ancient 
(~ 10 Gyr) initial episode of star formation. The only one 
exception to this result is IZw 18 (Izotov & Thuan l2004i 
located at ~ \5Mpc and believed to be relatively young 
(< 2 Gyr). In the appendix we briefly comment on this. 
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- We have discussed the spatial distribution of stars in dif- 
ferent age intervals. We find the distribution of young stars 
to be clumpy and associated with the regions of high emis- 
sion from the neutral gas, but not coincident with the max- 
ima of H i column density. Older main-sequence stars (with 
age between 0.7 and 2 Gyr) and intermediate-age (RC and 
RGB) stars show an extended, although not perfectly uni- 
form, distribution. There is no obvious correlation between 
the distribution of these relatively old stars and that of the 
neutral gas. 

- By considering the energy input from supernova explo- 
sions, we find that the "hole" in the Hi distribution cannot 
have been produced by a supernova driven wind unless, 1) 
the hole is unfeasibly old, 2) supernovae are able to im- 
part 2 orders of magnitude more kinetic energy to the ISM 
than previously thought, or 3) the IMF is very top-heavy. 
We can reasonably exclude 1 and 3, and 2 is very unlikely. 
By far the most difficult obstacle for the supernova wind 
driven hypothesis is the fact that the stellar surface bright- 
ness within the H i hole is lower than that outside the hole 
for any stellar age range. We find an appealing alternative 
to be that the H i ring (rather than "hole") is the result of a 
recent accretion of an intra-group H i cloud with a non-zero 
angular momentum. 
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Appendix A: On the age of I Zw 1 8 

The discovery of a population of presumably very old stars (red 
HB stars) in SagDIG argues against metal-poor galaxies being 
necessarily young, i.e. born in recent times. In this regard, we 
found it of interest to re-examine the properties of the prototype 
of candidate young galaxies, IZw 18. 

When comparing the two galaxies, it should be remem- 
bered that the metallicity difference between the two galax- 
ies is not large: IZw 18 is not so different after all from the 
most metal-poor Local Group dwarf galaxies like SagDIG or 
Leo A [see Aloisi et al. d2003t for a recent compilation of 
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Fig. 13. The mp8i4w, (»zf555w - twf814w) color-magnitude diagram of IZw 18 from ACS archival data. Panel A compares the 
photometry of Izotov & Thuan ( 120041 dark dots) with the result of our reanalysis (light dots), showing a better recovery of the 
fainter red stars. The HST/WFPC2 color-magnitude diagram of Sextans A (Dolphin et al. 2003) is shown for comparison in 
Panel B, together with hand-drawn fiducial lines marking the major evolutionary sequences (MS, blue and red super-giants, and 
RGB). A box delimits red giant stars in the RGB tip region. The same fiducial loci are superposed to the CMD of IZw 18 in 
Panel C, corrected for the distance and reddening of Sextans A and shifted to the distance of I Zw 18, assuming (m - M) a = 30.88 
(15 Mpc) with E(B - V) — 0.15. On the same assumptions, theoretical isochrones from Girardi et al. (2004, priv. comm.) with 
Z = 0.0004 and ages 7, 80, 281 Myr and 10 Gyr are over-plotted in Panel D. 



heavy elements abundances in IZw 18]. Yet IZw 18 is a key 
object in that it lacks a positive detection of a red giant branch, 
which seems to imply that the galaxy is younger than about 1 
Gyr. Even the recent HST/ACS deep study of Izotov & Thuan 
(2004 1 was unable to trace an "old" RGB. These authors place 
IZw 18 further out than previously thought, at 15 Mpc instead 



of 1 Mpc, and set an upper limit to the age of I Zw 1 8 between 
1-2 Gyr. 

By applying the reduction methods described in Sect. [5] to 
the ACS data of Izotov & Thuan (2004), we have been able to 
obtain a substantially improved CMD of IZw 18. The result of 
our reanalysis is presented in Fig. where we show a new 
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7«F8i4W> ( m F555w ~~ m F8i4w) CMD for a combined field contain- 
ing both the main body and the "C" component of I Zw 18. 

The new diagram is shown in Panel A, along with that of 
Izotov & Thuan (2004i. The comparison indicates an excellent 
match for the CMD features in common, with no appreciable 
photometric differences between the two reductions. Our re- 
duction however provides access to a higher completeness ev- 
erywhere, and deeper photometry. This result is probably to be 
attributed to the updated geometrical distortion correction and 
pixel oversampling technique employed here (see Anderson & 
King 2000) which appears to improve the effectiveness of PSF 
fitting. Note that the selection criteria on allstar parameters 
are similar for both reductions. 

It is for the redder population that the improvements in 
our diagram are most noticeable. In particular, we can iden- 
tify with some confidence the sequence of red stars having col- 
ors (mp555w ~ 7Wf814w) =1-2, made of red super-giants. This 
identification, crucial to set observational constraints on the re- 
construction of the SFH of IZw 18, is supported by compari- 
son with the CMDs of other dwarf irregular galaxies. Panel B 
shows the HST/WFPC2 color-magnitude diagram of Sextans A 
obtained by Dolphin et al. J2003t . where we have marked the 
main features (MS, blue and red super-giants, and RGB) with 
fiducial lines. The same fiducial loci are superimposed onto the 
CMD of IZw 18 in Panel C, assuming a distance of 15 Mpc for 
the galaxy. While the higher luminosity reached by MS stars in 
IZw 18 is consistent with its stronger star-burst, this compari- 
son shows an intriguing agreement for the older evolutionary 
sequences, in particular around the red super-giant and RGB 
tip regions. This conclusion is confirmed by a comparison with 
isochrones from Girardi et al. with Z = 0.0004 and ages 7, 
80, 281 Myr, and 10 Gyr, which are over-plotted on the data in 
Panel D. 

The termination point of the 10 Gyr isochrone overlaps 
with the clump of stars seen just fainter than rnpsnw ~ 27.3. 
Note also that the magnitude difference between the tip of the 
old RGB and that of the RSG is typical of dwarf irregulars, 
spanning 3.5^4.5 magnitudes in / with some dependence on 
the SFH of each galaxy. Although this result does not yet rep- 
resent a strong case for the identification of an "old" RGB, it 
seems very encouraging. Our reanalysis gives supports to the 
15 Mpc distance proposed by Izotov & Thuan (2004i, and pro- 
vides preliminary yet intriguing evidence for a RGB cutoff at 
'«F8i4w ~ 27.3. In this scenario, IZw 18 would not appear sub- 
stantially different from extremely metal-poor (yet old) Local 
Group dwarfs like SagDIG and Leo A. Efforts to push even 
deeper the capabilities of ACS would be worthwhile to confirm 
this hypothesis. 



